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Models of exchange-bias in thin films have been able to describe various aspects of this technologically 
relevant effect. Through appropriate choices of free parameters the modelled hysteresis loops adequately 
match experiment, and typical domain structures can be simulated. However, the use of these parameters, 
notably the coupling strength between the systems' ferromagnetic (F) and antiferromagnetic (AF) layers, 
obscures conclusions about their influence on the magnetization reversal processes. Here we develop a 2D 
phase-field model of the magnetization process in exchange-biased CoO/(Co/Pt) XJ1 that incorporates the 
10 nm-resolved measured local biasing characteristics of the antiferromagnet. Just three interrelated 
parameters set to measured physical quantities of the ferromagnet and the measured density of 
uncompensated spins thus suffice to match the experiment in microscopic and macroscopic detail. We use 
the model to study changes in bias and coercivity caused by different distributions of pinned uncompensated 
spins of the antiferromagnet, in application-relevant situations where domain wall motion dominates the 
ferromagnetic reversal. We show the excess coercivity can arise solely from inhomogeneity in the density of 
biasing- and anti-biasing pinned uncompensated spins in the antiferromagnet. Counter to conventional 
wisdom, irreversible processes in the latter are not essential. 

I n coupled ferromagnetic- ("F") and antiferromagnetic ("AF") thin films exchange-bias can arise for fixed AF 
I magnetic structures 1 . The effect, widely used in contemporary magnetic devices such as giant- and tunnel - 
I magnetoresistive thin-film sensors 2,3 , is set up when the F structure "imprints" a stabilizing structure in the AF 
upon cooling below the AF Neel temperature. It is manifested primarily as a lateral shift of the hysteresis loop of 
the F layer 4 . Often the width of the loop also increases with the onset of exchange-bias (e.g. 58 ), recently prompting 
studies of the use of this excess coercivity as a proxy for the degree of sub-monolayer Co-oxidation 9 . 

Microscopic model views of AFs and F-AF interfaces have provided insight into the mechanisms by which 
these features arise at the smallest scales 10 " 17 . It is clear that one characteristic of the materials' systems associated 
with exchange bias is the existence of pinned uncompensated spins antiparallel to the F-magnetization 18 " 21 for Co, 
Fe or permalloy coupled to CoO. The excess coercivity, on the other hand, has been circumscribed in the models 
to the effects of irreversible processes in the AF 16 . Understanding these phenomena at scales that reveal domain 
wall motion in polycrystalline films is important because many (if not most) applications rely on magnetization 
reversal through this process - or their impediment 21 " 24 . Hence the relevance of work by Fujiwara et al. and Stiles 
et al. 22,2526 that studied the influence of distributions of AF crystallite orientations and anisotropy in exchange bias. 
However, in these studies the evolution of domain walls in the F during reversal could not be accounted for 
because of the exceedingly high computational cost of modelling macroscopic systems. Other works have recently 
overcome this limitation 27 " 30 . But the role of the AF/F coupling, which is perhaps the least well understood 
component of exchange bias, is obscured by the reliance on free parameters for its description, and by the large 
variety of experimental results. 

With this work our goal is to avoid free parameters in a model description of F reversal processes in typical 
exchange bias systems. Instead, we want to rely exclusively on measured (or literature) values for the samples 
described, and show how accurate a model description can be, as assessed from measured microscopic domain 
images in applied fields and from hysteresis loops. Accordingly, the dispersion in the measured sample data used 
as model input must be small - hence we studied a single sample. In particular, on the contentious issue of the 
distribution of pinned uncompensated spins, we can use 10 nm-resolved experimental data from Schmid et al. 21 . 
As for the coupling, its average can be deduced from magnetometry of the exchange-bias field H ex . Various 
statistics of the anisotropy of the AF on a granular scale can be accessed with techniques such as that proposed by 
Vallejo Fernandez et al. 31 . This information would be essential for a correct account of temperature dependent- 
and training phenomena, and lacking it, we do not attempt to describe these effects (Incorporating them would be 
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possible, but beyond the scope of this work). The spatial distribution 
of F and AF anisotropy on the other hand cannot easily be furnished 
by experiment. So as not to shape the model outcome with our 
particular choice of distribution, we select the most general distri- 
bution possible, a Gaussian one. Improved models would base the 
anisotropy distributions on statistics of measured Barkhausen ava- 
lanches in the sample 32 , or infer spatial distributions from the mar- 
ginal changes of domain boundaries with applied field at room 
temperature, where there is no exchange bias. 

Our model is a 2D phase-field model, similar to those used to 
describe ferromagnetic films in relation to the role of disorder in 
domain dynamics. For instance, return point memory effects 33 , 
Barkhausen avalanche distributions 32 , and the role of defects in the 
domain reorientation under the influence of an oscillating external 
field 34 have been investigated in this fashion. These studies did not 
attempt to match experiment quantitatively, in part because the local 
domain pinning strength was not known; nor have they been imple- 
mented in the context of exchange-bias. From this perspective, our 
model differs from conventional ones in three important ways: First, 
it describes an exchange-bias system on a scale relevant for magnet- 
ization reversal processes governed by domain wall motion. Second, 
it incorporates an experimentally determined 21 10 nm-resolved dis- 
tribution of F-AF pinned uncompensated spins ( pin UCS) over the 
2 fim X 2 fim area modelled. As we will show in the following, the 
inhomogeneity of pinned UCS affects H ex and H c 13 ' 21,26 . Hence 
the importance of using experimental values for pinned UCS, and 
the implied refinement over previous models. And third, the model 
agrees quantitatively with the experimentally determined hysteresis 
loop and 10 nm-resolved domain evolution, using for both scales 
one and the same set of material parameters that are in agreement 
with commonly accepted experimental values. 

Results 

Model construction. Figure 1 illustrates schematically our CoOi nm / 
Co 0 .6 nm[Pt 0 . 7 nm/Coo.4 nm ]x2o exchange-biased system and the 
elements of the model counterpart. In this work, we obtain the 
domain dynamics from the Landau-Lifshitz-Gilbert equation 
(LLG) 35-37 governing the damped precession of the magnetization 
M of the ferromagnetic layer in presence of a field B: 

cm y 



mx [B + £(mxB) 



m 



where y is the gyromagnetic ratio, £ the Gilbert damping para- 
meter, and m(R, t) = M(R)/M S . M s is the saturation magneti- 
zation, assumed uniform. B is the magnetic field. We write 
B= — l/M 5 (57Y[m]/c)m + Q(R, f), using the functional derivative 
of the hamiltonian H and a gaussian stochastic process Q(R, t) 
which accounts for the finite temperature effects. ((Q(R, t)) = 0 
and <Q(R, f)Q(R', f')> = <$(f - t')S(R - R')2k B W M s y; k B is the 
Boltzmann constant and Tthe temperature.). 

We assume uniform and small thickness d and perpendicular 
magnetization. The latter is represented by a scalar dimensionless 
field m(x, y, t) such that m(r, t) = m(x,y, t)z. This assumption 
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Figure 1 | Schematic of our exchange-biased CoOx nm /Co 0-6 nm [Pto.7 nm / 
Co 0 .4 nml x2o multilayer. 



implies the description will be accurate as long as the domain wall 
is much narrower than the domain width 37 " 39 . 
The system's hamiltonian can then be written as: 
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The first term in (2) represents the anisotropy energy given the 
uniaxial anisotropy K U (R). 

The second term in Tt represents the exchange interaction in the F 
layer, described by the exchange stiffness A. 

The third term represents the long-range, non-local, stray field (fi 0 
is the vacuum permeability) for small film thickness. 

The fourth term contains two fields. On the one hand there is an 
external uniform field H ext whose value is quasi-statically ramped up 
and down in time; On the other, an effective field H ucs (x, y) = 
H e bP p (x, y)/(pp), describing the local biasing effect in terms of the 
macroscopic, measured exchange-bias field H eb and the measured 21 
pin UCS density p p (x,y) (cf. Fig. 2(g)). Simulations of the field-cooled 
hysteresis loops require using H s ^ s (x, y) as described in Fig. 2(h) 
(see See Supplementary Information for details on its construction 
from H ucs (x, y)). 

Substituting (2) in the over-damped limit £>4 of eq. (1) and in the 
approximation of thin domain walls 38 , the equation for the F domain 
dynamics becomes: 
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Figure 2 | Model validation, (a) Experimental 10 K domain pattern at 
0 mT applied field (from Ref. 21). (b) and (c) Corresponding domain 
patterns at 100 mT and 200 mT applied field, (d), (e) and (f) Simulated 
domain patterns, (g) Experimental pattern /z^es usec ^ f° r (d)-(e). (h) 
Modified h ucs pattern (h s ^ cs ) for the simulation of hysteresis loops, 
obtained from (g) through inversion of the areas corresponding to the light 
domains of (d). (i) Comparison of hysteresis loops from experiment 21 with 
the model result. 
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(3) 

+ ^W-^r/cs(r) + ^(r, T)J+^m, 

where we have introduced explicitly the dimensionless units r = R/d, 
t = ty£fi 0 M s , h ext = JWM 5 , /z^cs = H ucs /M s and g(r, t) = Q(R, f)/ 
^ 0 M 5 . The dimensionless constants oc= (K u ) / /n 0 M 2 and 
P = A/in 0 M^d 2 are the reduced anisotropy and exchange stiffness, 
respectively. p(r) is the distribution that accounts for the anisotropy 
inhomogeneity, and is not exactly known. To avoid introducing 
overly restrictive assumptions into the model we choose p(r) to be 
Gaussian-distributed spatially-uncorrelated noise with (p) = 0 and 
variance r], i.e. (p(r)p(r r )) = rj3(r — r'). Together, a, ft and rj consti- 
tute the only parameters of the model. 

Model validation with experiment. In Fig. 2 we compare the model 
results for a = 6.6, /? = 0.14 and rj = 1.88 X 10" 4 with experimental 
results from 21 . Figures 2 (a)-(c) are the experimental, 10 K domain 
structures of the zero-field cooled system for 0 mT, 100 mT and 
200 mT applied field. At 0 mT Fig. 2 (d) is the stable domain 
structure calculated with the model starting with the domain 
structure from 2 (a). The detailed resemblance between Fig. 2 (a) 
and (d) shows the measured pattern is stable for the model, given the 
magnetization dynamics it describes. Starting from Fig. 2 (d) and 
gradually increasing the applied field, and letting the system evolve 
using Eq. 3 to a stable domain structure results in Figs. 2 (e) and (f). 
The shape of the simulated domains matches the experimental one, 
and a high level of detail is reproduced by our model, although some 
discrepancies are apparent. 

We then simulate a magnetometry measurement at 10 K after 
cooling in a 1 T field, that is, a typical exchange-bias measurement, 
Fig. 2 (i). This we do using exactly the same values of a, and rj 
employed for the above domain evolution. However, we cannot use 
exactly the same distribution of pinned UCS as before, since it is set 
by the F's magnetization structure during cooling, which is now 
different. Neither can we directly measure the pinned UCS from this 
experiment preparation 40 : Recall that a uniformly magnetized F film 
would produce no stray field for the MFM to detect, so in the field 
cooled case the MFM would only image the inhomogeneity of the 
pinned UCS on a local scale, but not their average density. But note 
that the F, over the areas inside its domains, sets the pinned UCS in 
the same way as an F saturated with the appropriate orientation 
would. Thus for the loop simulation we construct the map of pinned 
UCS shown in Fig. (h) from the one obtained after zero-field cooling, 
shown in Fig. (g) by inverting the pUCS under the white ferromag- 
netic domains. 

The important loop characteristics agree well with experiment. 
Specifically, we observe a mean coercivity of 86.20 mT and an 
exchange field of 12.35 mT in the field- cooled case, which compare 
well with experiment (77.7 mT and 13.6 mT, respectively). 
Prominent features of the magnetization loop, such as the knee at 
the nucleation field and the subsequent more protracted approach to 
saturation are also displayed by the model results for the hysteresis. 
In the zero field cooled case, where we expect no macroscopic 
exchange bias, we obtain a match of comparable quality as for the 
field-cooled case. For the corresponding comparison at 300 K, also in 
Fig. 2 (i), we use a = 5.8 (/? and rj retain their 10 K-values), thus 
accounting in a qualitative manner for the reduced anisotropy at 
higher temperatures. Furthermore, consistent with the absence of 
*"'"UCS above the AF blocking temperature, we set h UC s = 0 in this 
case. As expected, the simulated loops are symmetric and have a 
considerably reduced coercive fields. 

Contributions to domain dynamics from anisotropy- and coupling- 
inhomogeneity. Our model allows us to investigate how the domains 



evolve greater detail, for example looking at intermediate field levels 
(See Supplementary Information for a video of the simulated F 
domain evolution and detailed images of domain patterns.). Also, 
for instance, we can investigate the changes in domain evolution in 
the hypothetical case of h UC s = 0, that is, when domain boundary 
pinning is controlled solely by anisotropy inhomogeneity in the F 
layer and there is no net exchange- coupling between F and AF. 
Figures 3(a)-(c) show the resulting domain patterns at 0, 100, and 
200 mT in black and white. For comparison, we superimpose in 
yellow the contours of the corresponding domain boundaries from 
experiment. Likewise, we study the converse case of hues = ^ucs anc ^ 
rj = 0, i.e. K U (R) = (K u ), when pinning from anisotropy 
inhomogeneity is negligible. We show the results in Figures 3(d)- 
(f) (cf. Supplementary Information for a side-by- side summary of the 
images). 

Exchange-bias and coercivities for different inhomogeneity in 
local density of pinned uncompensated spins. With our model it 
is also straightforward to calculate the magnetization process 
characteristics for hypothetical distributions of pinned uncompen- 
sated spins hues, i.e. varying degrees of F-AF coupling inhomo- 
geneity. In Tab. 1 we summarize the coercivity and exchange-bias 
fields from experiment and the full simulation, as well as values 
reported by simulations where h ucs has been modified arbitrarily. 
The first two reference data sets comprise the simulation with 
hues = ^ucs discussed in Fig. 2 (See also Supplementary Informa- 
tion), and a further simulation setting hucs = ®> which confirms 
that the resulting loop is centered around 0 mT. It has a slightly 
reduced coercivity, and there is no exchange-bias in this case: 
Kb = 0. 

Next, we can set h U cs = h s u lt cs + (h s ^ cs ) to simulate a film with 
double the average density of h s ^ cs and with the same inhomogeneity. 
In other words, a film with greater average pin UCS but comparable 
levels of coupling frustration and local variations in the pinning 
ability of the AF. The resulting hysteresis curve is laterally shifted 
by an amount twice as large as H eb in the first simulation, confirming 
the experimental results based on microscopic domain pattern 
evolution. Notably the enlarged coercivity, which typically accom- 
panies exchange-bias, remains at the original levels. 

These findings change if, conversely, we set hues = 2 x h s ^ cs — 
(hues) t0 simulate a film with the original average density of 
pin UCS but twice the amplitude of the inhomogeneity. In that case 




Figure 3 | Simulated 10 K magnetization patterns (black & white) for 
different hypothetical model inputs, put in comparison with experiment 
(yellow trace), (a) Simulated domains at 0 mT applied field using rj = 1.88 
X 10" 4 and h ucs = 0. (b) Idem at 100 mT. (c) Ibidem at 200 mT. (d) 
Simulated domains at 0 mT applied field using rj = 0 and hues = h^s from 
Fig. 2(g). (e) Idem at 100 mT. (f) Ibidem at 200 mT. In this simulation the 
sample already saturates at this field level. 
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the model results in a greater coercivity without significant changes 
in the exchange-bias. 

Discussion 

From Figure 3 we can see that simulations without anisotropy 
inhomogeneity or without pinned uncompensated AF spins cannot 
match the experimental domains as accurately as the full simulation 
with h ucs = h e ^ s and n = 1.88 X 10" 4 , Figs. 2(d)-(f). Despite such 
deficiencies, the simulation with h UC s = 0 tracks the experiment with 
reasonable accuracy for 0 and 100 mT applied field, and shows more 
prominent deviations from it only at 200 mT. This high field beha- 
vior is compatible with the smaller number of energy minima of 
sufficient strength in the absence of coupling to the AF. 

As a means to control domain boundary motion h ucs appear to be 
slightly less effective than K u inhomogeneity, which agrees with 
H ex < H c . In particular, Figs. 3(d)-(f) depart markedly from the other 
simulations already at 0 mT applied field. Furthermore at 200 mT 
the simulation would predict magnetization saturation, so, clearly, 
pinning from h ucs alone is unable to describe the measured domain 
structure at 200 mT. These observations do not imply a subordinate 
role of p/n UCS in exchange-bias; on the contrary. They show that even 
without anisotropy inhomogeneity to pin the F domain walls a mag- 
netization structure is retained up until at least 100 mT, which would 
not be possible if p, '"UCS did not pin the domains. 

Moreover, the discrepancies between model and experiment 
depend significantly on the local values K U (R), the distribution of 
which we have carefully kept in the most generic form. Because of 
that, K U (R) very likely differs from the real distribution of K w and 
simulation inaccuracies are to be expected. Nevertheless, the model 
captures the essential physics of the magnetization reversal in the 
presence of domain wall motion. An anisotropy distribution with a 
sufficiently large number of free parameters could of course be 
adjusted so as to yield a more perfect match between model and 
experiment. However, this strategy would obscure, rather than clarify 
the magnetization reversal mechanism here. 

We can gain additional insight into how the reversal is affected by 
the pz "UCS and K u distributions when we look at simulations of the 
magnetization loops. For instance, we can reexamine the link 
between exchange-bias and average pin UCS found in microscopic 
observations 21 . Hysteresis loops at 10 K are the macroscopic 
counterpart to these observations, which could not have been carried 
out in experiment due to the impossibility of arbitrarily changing the 
density of pinned uncompensated spins. 

Table 1 confirms that exchange-bias is roughly proportional to the 
average density of pinned uncompensated spins. In particular, 
pinned uncompensated AF spins aligned parallel to the F are det- 
rimental to exchange bias. Moreover, the inhomogeneity of the 
pin UCS leads to an excess coercivity. The importance of this finding 
is that rotating UCS and irreversible processes in the AF are not a 
necessary condition for excess coercivity, as the prevalent thinking 
holds 16 . Nor is their role in F reversal at low temperature expected to 
be major, except perhaps in the case of strongly coupled 41 " 43 F-AF 
layers. 



The picture that emerges is one where the inhomogeneity of the F's 
anisotropy and of the AF's pinned uncompensated spins largely 
determines the details of the F reversal. From our data 21,23,44 ' 45 as well 
as from XMCD measurements 19 it is clear that the pinned uncom- 
pensated spins exist in exchange bias systems in relatively high areal 
densities, exceeding —10% of a monolayer. This is in apparent con- 
tradiction to the observed exchange fields, which are far smaller than 
would be expected if all pinned UCS coupled with bulk- order 
exchange constants to the F (note that by measuring pinned UCS 
we consider only the part that is stable in applied fields). The problem 
stems from the default assumption that all pinned UCS participate in 
the coupling to the F. Notice though, that the pinned UCS found in 
the aforementioned experiments align antiparallel to the F moments. 
Consequently they cannot be aligned directly by the cooling field. 
Instead, these pinned UCS must be exchange-coupled to the F (anti- 
ferromagnetically, possibly via a superexchange mechanism), or 
otherwise exchange- coupled to such UCS. In either case it follows 
that these pinned UCS must be located at or near the F-AF interface. 
But of them, only the ones that couple directly (antiferromagneti- 
cally, as discussed) to the F generate an exchange bias effect. The 
remaining ones do not provide additional coupling, explaining the 
weakness of the exchange bias in the presence of a surprisingly high 
density of pinned UCS. Our conclusion is supported by reflectometry 
experiments that revealed the pinned UCS existed over a film thick- 
ness larger than the roughness of the interface 19 . Hence models 
assuming a sharp interface or uncompensated spins located solely 
at an atomically sharp interface between an F and AF seem inappro- 
priate to explain exchange bias. 

Thus our 2D phase-field model of exchange-bias systems based on 
general assumptions and parameters set by experiment, matches 
experiment in macroscopic and microscopic detail, and on this basis 
is able to establish that: 1) The magnitude of the average p/ "UCS 
density determines the exchange -bias field H eb in spite of the fact 
that only a part of these p/n UCS couple directly to the F. 2) The spatial 
inhomogeneity of the p/n UCS governs the evolution of the domain 
pattern on a local scale and gives rise to excess coercivity associated 
with exchange bias. 3) Irreversible AF processes and^"UCS rotation 
need not be invoked to explain excess coercivity. 4) The average 
coupling between the F moments and the pin UCS is weak compared 
to intrinsic coupling constants describing the exchange in ferromag- 
nets and antiferromagnets. 

Methods 

It is important to discuss the extent of the adjustments admitted in the course of 
matching the model results to experiment (See Supplementary Information for a table 
of material parameters and adjustment guidelines). We emphasize that a, P and rj are 
not arbitrary ("free"). Specifically, M s and d are precisely known for the particular 
film we used 21 . With regard to a, because (K u ) is not available for our film the model 
calculations are carried out using a literature value 46 subject to fine adjustments of few 
%. This is acceptable given K u may differ slightly for two films of the same nominal 
characteristics and fabrication process. We restrict the RMS inhomogeneity of K u to 
not more than about 20%, effectively limiting the range of possible values of rj to small 
positive numbers of order 10~ 4 . Determining P further requires specifying A. We use 
bulk values found in literature 46 and apply corrections for the different dimensionality 
of our model, as called for by data on the domain wall width 3 dw in our system 46 and 
the relation 3d w = ndy/ 'P/oc for Bloch- walls. The values oc ini = 6.25, p ini = 0.135 and 
rj ini = 1.5 X 10" 4 , used for the first model calculation on the basis of which further 
adjustments follow, are determined in this way. 

The simulation of the domain pattern evolution over a series of applied fields 
consists in integrating Eq. (3) numerically (in Fourier space, to circumvent the non- 
locality. Cf. Supplementary Information for technical details) starting from the 
known experimental 21 zero -applied- field F-domain structure (at T = 10 K), Fig. 2(a). 
Hysteresis loops follow trivially from the series of patterns. Taking into account the 
influence of a, P and rj on the domain evolution, we modify them slightly so that the 
zero applied field magnetization pattern in Fig. 2(a) becomes a stationary state of our 
model. Even smaller manual adjustments of a, P and rj follow, to ensure the simulated 
evolution of the domain pattern with the applied field about matches the other 
available experimental domain patterns, i.e. Fig. 2(b) -(c), at h ext =100 and 200 mT, 
respectively. A final manual fine-tuning of parameters seeks a fit of the experimental 
10 K hysteresis loop, Fig. 2(i). Note that an automated search for best fitting {a, p, rj] 
values is possible (e.g. relying on crosscorrelations for domain patterns as a fitness 
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function 33 ) but not practical at the moment. It would not change the model conclu- 
sions substantially. 

By this process we arrive at a = 6.6, /? = 0.14 and n = 1.88 X 10~ 4 , which depart 
only slightly from the initial values. 
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